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Analytic  and  simulation  results  are  presented  which  describe  the 
statistics  of  the  duration  and  separation  of  scintillation  which  results 
when  radio  frequency  signals  propagate  through  randomly  ionized  media. 
This  work  is  applicable  to  the  problems  of  satellite  communication  and 
space  based  radar  observation  through  disturbed  ionospheric  channels  that 
result  from  high  altitude  chemical  release  or  nuclear  detonation.  In 
these  environments,  a  radio  frequency  signal  that  traverses  the  disturbed 
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20.  ABSTRACT  (Continued) 

Ionosphere  a  single  time  has  Rayleigh  amplitude  statistics  and  power 
spectral  densities  that  range  from  f-3  in  the  slow  fadinq  limit  to 
Gaussian  in  the  fast  fading  limit. 

The  amplitude  statistics  of  the  received  signal  are  extended  to  mono¬ 
static  and  bistatic  propagation  through  the  disturbed  ionosphere.  The 
bistatic  geometry  is  applicable  to  either  a  bistatic  radar  or  to  a  two-way 
communication  link  that  utilizes  a  linear  transponder.  The  statistics  of 
the  amplitude  used  for  demodulation  in  dual  channel  communication  links 
are  presented  for  systems  that  use  spatial,  frequency,  or  other  diversity 
techniques.  The  distributions  of  the  received  amplitude  are  used  to 
calculate  bit  error  rates  for  two-way  or  dual  charnel  communication  links 
and  for  various  PSK  and  FSK  demodulation  techniques.  These  results  are 
compared  with  bit  error  rates  in  single  channel  or  one-way  communication 
links  and  with  bit  error  rates  in  communicaL  In  links  through  the  ambient 
ionosphere. 

Analytic  expressions  are  obtained  for  the  mean  number  of  level  cross¬ 
ings  and  the  mean  duration  and  separation  of  scintillation  above  or  below 
arbitrary  power  levels  in  the  one-way,  monostatic,  and  bistatic  propaga¬ 
tion  geometries.  These  mean  quantities  are  plotted  for  a  Gaussian  power 
spectral  density  or,  equivalently,  for  a  fast  fading  environment.  The 
mean  duration,  separation,  and  number  of  level  crossings  for  arbitrary 
fT  power  spectral  densities  can  be  obtained  from  the  Gaussian  values  by 
scaling  factors  which  are  calculated  for  jj  in  the  range  3  to  10. 

The  distributions  of  the  number  of  level  crossings  and  the  duration 
and  separation  of  fast  fading  scintillation  are  obtained  using  Monte  Carlo 
techniques  for  the  one-way,  monostatic,  and  bistatic  propagation  geome¬ 
tries  and  for  dual  channel  communication  links.  The  simulation  results 
are  presented  in  tabular  form  and  are  compared  with  analytic  results  where 
possible.  Histograms  of  the  duration  and  separation  of  scintillation  for 
the  one-way  propagation  geometry  are  also  presented.  Qualitative  compari¬ 
sons  are  made  of  the  effects  of  Rayleigh  fading  conditions  on  ore-way 
versus  two-way  and  single  channel  versus  dual  channel  communication  links 
or  monostatic  versus  bistatic  radar  systems. 
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SECTION  1 
INTRODUCTION 


The  performance  of  a  space  based  communication  or  radar  system 
operatinq  in  a  nuclear  environment  Is  deqraded  by  the  result inq  random 
amplitude  and  phase  fluctuations  of  the  received  radio  frequency  slqnal. 

In  qeneral,  the  random  amplitude  fluctuations  produce  an  Instantaneous 
slqnal  power  that  is  less  than  the  averaqe  power,  while  occasionally,  the 
Instantaneous  slqnal  power  will  be  larqer  than  the  averaqe  power.  The 
first  situation  Is  referred  to  in  this  report  as  a  fade  and  the  latter  as 
a  flare.  Because  the  propaqation  disturbances  are  correlated  In  time, 
these  fades  or  flares  have  a  finite  mean  duration  and  mean  separation. 

More  qeneral ly,  fades  below  an  arbitrary  level  or  flares  above  that  level 
have  varyinq  probability  of  occurence,  mean  duration,  and  mean  separation 
as  the  level  is  chanqed  relative  to  the  mean  received  power.  It  Is  the 
purpose  of  this  report  to  describe  the  statistics  of  the  duration  and 
separation  of  fades  and  flares  for  several  cases  of  Importance  to  space 
based  communication  and  radar  systems.  This  work  was  performed  In 
response  to  requests  from  contractor  and  SPO  personnel  Involved  In  the  DSP 
and  IONDS  proqrams. 

Three  propaqation  qeometrles  are  considered  and  are  denoted  one¬ 
way,  monostatic,  and  bistatic  throuqhout  this  discussion.  The  one-way 
case  corresponds  to  a  one-way  communication  link  and  the  monostatic  case 
corresponds  to  a  monostatic  space  based  radar.  The  bistatic  case  corre¬ 
sponds  to  either  a  bistatic  space  based  radar  or  to  r  two-way  communica¬ 
tion  link.  The  bistatic  case  reduces  to  the  one-way  case  In  the  event 


that  only  one  of  the  propagation  paths  Is  throuqh  a  Rayleigh  fadit  q  chan¬ 
nel  while  the  other  path  is  through  an  undisturbed  channel.  However,  it 
is  expected  that  this  situation  has  a  low  probability  of  occur ance  and  it 
is  mentioned  here  for  completeness  only.  Dual  channel  communication 
systems  which  combine  the  received  siqnal  from  two  independent  one-way 
communication  channels  before  demodulation  are  also  discussed  in  this 
report.  In  practice  the  independent  channels  can  be  achieved  by  the  use 
of  spatial  diversity  (e.q.,  separated  receiver  antennas  or  transmitter 
locations)  or  other  diversity  techniques  such  as  multiple  transmitter 
frequencies.  Two  different  algorithms  for  combining  the  siqnals  into  the 
amplitude  used  for  demodulation  have  be&n  considered.  The  summation  algo¬ 
rithm  adds  the  siqnal  amplitudes  from  the  two  channels  while  the  maximum 
alqorithm  selects  the  larqer  of  the  two  amplitudes. 

The  first  and  second  order  statistics  pf  the  received  signal  for 
the  monostatic  and  bistatic  propagation  cases  arp  described  in  the  second 
section  assuminq  Rayleigh  amplitude  statistics  and  Gaussian  temporal  auto¬ 
correlation  functions  or  equivalently  Gaussian  ptjiwer  spectral  densities 
(PSD)  for  the  received  one-way  voltaqe.  The  distribution  of  the  amplitude 
used  for  demodulation  is  also  calculated  in  this  section  for  the  two  dual 
channel  communication  link  algorithms  and  is  compared  with  a  sinqle  chan- 

I 

nel  or  one-way  communication  link.  As  an  application  of  the  first  order 
statistics  results,  binary  error  probabi  lit ies  are  calculated  in  the  third 
section  for  two-way  and  dual  channel  communication  links  under  slow 
Rayleiqh  fadinq  conditions  and  for  various  frequency-shift  keying  and 
phase-shift  keyinq  demodulation  techniques. 

The  mean  number  of  level  crossings  per  unit  time  for  Rayleiqh 
fadinq  with  an  arbitrary  PSD  is  calculated  next  in  closed  form  for  the 
three  propagation  cases.  The  mean  duration  and  separation  of  fades  and 
flares  is  obtained  from  the  mean  number  of  level  crossings.  These  anal¬ 
ytic  results  provide  a  qualitative  comparison  of  the  effects  of  Rayleiqh 

10 
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fading  conditions  on  one-way  versus  two-way  communication  links  or  mono¬ 
static  versus  bistatic  radar  systems.  The  mean  number  of  level  crossings 
and  the  mean  duration  and  separation  of  fades  and  flares  are  plotted  for  a 
Gaussian  PSD.  These  numerical  results  can  be  scaled  to  obtain  the  mean 
Quantities  for  an  arbitrary  f*M  PSD  using  scaling  factors  which  are 
calculated  in  Appendix  A. 

The  problem  of  calculating  the  probability  density  functions  of 
tiie  duration  and  separation  of  fades  or  flares  cannot  be  solved  in  closed 
form  for  even  the  simplest  propagation  geometry.  Thus  Section  5  briefly 
outlines  the  method  of  generating  random  seauences  of  voltages  with  the 
desired  first  and  second  order  statistics.  The  distributions  of  the  dura¬ 
tion  and  separation  of  fades  and  flares  for  a  Gaussian  PSD  are  then 
obtained  in  numerical  form  from  the  random  data.  These  simulation  results 
are  tabulated  in  Section  6  and  are  compared  with  the  analytic  mean  quanti¬ 
ties.  The  one-way  results  are  compared  in  Appendix  C  with  mean  fade 
duration  and  separation  data  from  several  multiple  phase  screen  (MPS) 
real izat ions .  Histograms  and  cumulative  distributions  of  fade  and  flare 
separation  and  duration  are  presented  in  Appendix  D  for  the  one-way  prepa¬ 
ration  case. 

The  major  conclusions  of  this  study  are  summarized  in  the  final 
section.  They  are: 

1)  To  achieve  a  demodulator  bit  error  rate  of  10*3  under  slow 
Rayleigh  fading  conditions,  a  two-way  communication  link*  requires  about  9 
dB  more  bit  energy-to-no ise  density  than  does  a  one-way  link  while  a  dual 
channel  communication  link  requires  about  10  dB  less  bit  energy-to-noise 
density  than  does  a  single  channel  link. 


*  Assuming  no  satellite  processing  and  a  perfect  (linear)  transponder. 
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2)  The  mean  duration  of  fades  below  the  mean  power  is  longest 
for  the  monostatic  geometry  and  shortest  for  the  one-way  geometry  with  the 
bistatic  geometry  mean  fade  duration  falling  between  the  other  two.  For 
fades  more  than  11  dB  below  the  mean  power,  the  monostatic  case  also  has 
the  shortest  mean  separation  between  fades. 

3)  When  compared  to  fades  in  single  channel  communication  links 
which  are  more  than  10  dB  below  the  mean  power,  the  10  dB  fades  of  dual 
channel  links  are,  on  the  average,  two  thirds  as  long  and  spaced  three 
times  farther  apart. 


SECTION  2 

FIRST  AND  SECOND  ORDER  STATISTICS  OT  THE  RECEIVED  SIGNAL 


In  this  section  well  known  empirical  descriptions  are  presented 
of  the  first  and  second  order  statistics  of  the  one-way  received  voltaqe 
in  the  stronq  scatter inq  limit.  From  these  results,  the  first  and  second 
order  statistics  of  the  received  voltaqe  for  the  monostatic  and  bistatic 
cases  are  constructed.  The  distribution  of  the  received  power  for  the 
three  propaqation  cases  are  then  calculated  as  are  the  distributions  of 
the  combined  power  for  the  two  dual  channel  algorithms. 

2.1  STATISTICS  OF  THE  ONE-WAY  VOLTAGE 

After  propaqatinq  one-way  throuqh  an  ionized  layer,  the  received 
voltaqe  can  be  expressed  as  two  independent  quadrature  components 

P(t)  =  x(t )  ♦  iy(t)  (2-1) 

where  p(t)  is  the  received  complex  voltaqe  with  in-phase  component  x(t) 
and  quadrature-phase  component  y(t).  Both  MPS  calculations1.2  and 
theoretical  calculations  in  the  stronq  scattering  reqime3  show  that  trans¬ 
mission  of  an  initially  constant  amplitude  and  phase  siqnal  one-way 
throuqh  a  stronqly  turbulent  layer  results  in  a  received  siqnal  whose 
quadrature  components  x(t)  and  y(t)  are  independent  Gaussian  random 
variables  with  joint  probability  density  function 

f(x,y)  =  exp[-(xz+y2)/2cr2]/2Tto2.  (2-2) 

(--  <  x  <  <  y  <  -  ) 
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At  two  different  times,  the  in-phase  components  of  p( tx )  and  p(t2), 
denoted  X!  and  x2  respectively,  have  the  joint  Gaussian  probability 
density  function 

f ( x ! ,x2 )  *  exp[ -(xj -2pxt x2+x2 )/2o2 ( 1-p2 ) ]/2*oz/l-p2  .  (2-3) 

The  joint  probability  density  function  of  the  quadrature-phase  components 
yj  and  y2  is  identical  to  Equation  2-3.  The  autocorrelat ion  of  the  com¬ 
plex  voltage  is  then 

<p( t)p*( t+r )>  =  <x( t )x( t+x )>  +  <y( t )y( t+T )>  =  2o2p(t).  (2-4) 
In  the  strong  scattering  limit 

p (t )  =  exp[ -(t/tj )2)]  (2-5) 

where  Tl  is  the  one-way  signal  decorrelation  time.  General  f-M  power 
spectral  densities  with  corresponding  Bessel  function  forms  for  p(-r)  are 
considered  in  Appendix  A.  In  subsequent  developments,  the  strong  scatter¬ 
ing  limit  Gaussian  form  of  p(x)  will  be  used  as  an  example  with  appropri¬ 
ate  scaling  factors  included  to  render  the  results  valid  for  any  properly 
behaved  PSD.  It  will  be  seen  later  in  this  section  that  only  in  the  case 
of  one-way  propagation  through  the  ionized  medium  does  the  decorrelation 
time  of  the  received  signal,  t0,  equal  the  one  way  decorrelation  time 

*i  • 


The  amplitude  r  of  p(t), 
r  =  /x2+y2 


(2-6) 


has  a  Rayleigh  probability  density  function 


I 


and  the  average  power  of  the  signal  is 


<p(t)p*(t)>  =  <r*>  =  2oz  . 


(2-8) 


2.2 


GENERAL  FIRST  ORDER  STATISTICS 


The  received  voltage  q(t)  of  a  radio  frequency  signal  which 
traverses  an  ionized  layer  either  once  or  twice  can  be  represented  as 


q(t)  = 


p(t)  one-way 

p2(t)  monostatic 

Pj  (t)p, (t)  bistat ic 


(2-9) 


when  it  is  assumed  that  in  the  radar  cases  the  target  cross  section  is 
non-fluctuating  and  that  in  the  two-way  communication  link  case  there  is 
no  on  board  satellite  processing  and  the  transponder  used  is  perfectly 

i 

linear.  The  effect  of  target  cress  section  fluctuations  on  the  received 
signal  has  been  discussed  in  detail  elsewhere.1*  In  the  bistatic  case, 

Px ( t )  represents  the  signal  transmitted  along  path  1  and  incident  on  the 
transponder  or  the  radar  target  and  p2(t)  represents  the  signal  retrans¬ 
mitted  or  reflected  by  the  target  and  received  along  path  2.  The  received 
power  S  is  then  given  by 


S  =  q( t)q*( t)  = 


n  one-way 

r1*  monostatic 

2  2 

r1r2  bistatic 


(2-10) 


The  probability  density  function  of  S  is  easily  calculated  from  the 
Rayleigh  probability  density  function  of  r  as 


exp[ -S/<5>] /<S> 


f(S)  =  {  exp[  -V 2S/<S>] // 2S< S> 


2Kq[  2/  S/< S> ]  /< S> 
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one-way 
monostat ic 
bistatic 


(2-11) 


'I  ' 

\  , 

\  - 


t 


where  K0  is  a  modified  Bessel  function.  The  mean  power  is 


( 2°2 

one-way 

<s>  = 

!  8o4 

monostatic 

bistatic 

(2-12) 

where  o2  is  the  variance  of  the  underlying  Gaussian  probability  density 
function  of  Equation  2-2. 


The  cumulative  distribution  function  of  the  received  power, 

defined  as 

S 

F(S)  -  Prob  (received  power  _<  S}  s  /  f(S')dS‘  ,  (2-13) 


gives  an  indication  of  the  probability  of  occurrence  of  deep  fades  which 
are  detrimental  to  the  performance  of  communication  or  radar  systems. 
This  function  is  readily  calculated  for  the  three  cases  considered  with 

the  results 


-  exp[-S/<S>] 


('  _ 

F(S)  *  j  1  -  exp[-/2S/<S>] 

(  1  -  2/s7<S>  Kl[2/S/<S>]  bistatic 


one-way 

monostatic 


where  Kx  is  a  modified  Bessel  function. 

The  median  power  M,  defined  by 
F(M)  =  1/2, 


is  equal  to 


tn(2)  =  -1.6oB  one-way 

M/<S>  *  {  [tn(2)]z/2  =  -6.2dB  monostatic 

u  /4  =  -4 .OdB  bistatic 

o 


(2-14) 


(2-15) 


(2-16) 
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the  highest  probability  of  deep  fades,  the  one-way  case  has  the  lowest 
probability,  and  the  bistatic  case  falls  between  the  other  two.  For 
example,  the  one-way  case  has  a  probability  of  lO'*4  of  fades  of  40  dB  or 
greater  while  this  probability  is  about  8  *  10~4  for  the  bistatic  case  and 
is  about  1.5  x  10~2  for  the  monostatic  case. 


2.3 


GENERAL  SECOND  ORDER  STATISTICS 


The  frequency  of  occurrence,  the  duration  and  the  separation  of 
fades  or  flares  is  determined  by  both  the  first  order  or  amplitude  statis¬ 
tics  described  above  and  by  the  second  order  or  correlation  properties  of 
the  received  voltage  defined  in  Equation  2-4.  The  correlation  function 
for  the  monostatic  case  is 

<q( t)q*( t+x )>  =  <p2(t)p2*(t+x)> 


*  <*i>4  +  *1*2  +  4xl *2*1*2  “  *1*2  ’  *2*1  *  (?-19) 


+  2i<xlx2y1  -  Xj yt  y2  -  XiX2.y2  +  x2yfy2>  . 


Noting  that  x  and  y  are  independent,  zero  mean  random  variables,  the  imag¬ 
inary  part  of  <q(t)q*(t+x  )>  is  identically  zero.  In  addition,  it  is  easy 
to  demonstrate  the  following  using  Equation  2-3: 


.  2  2  s 
<Xj  x2  > 


<y2i  y\  >  *  a4  ( l+2p2  ) 


<*l*2*l*2>  =  <*l*2x*l*2>  =  (C2P)2 


(2-20) 


<xi  yf  >  =  <x2  *i  >  =  <**></>  = 


2  2. 


so  that  for  the  monostatic  case. 


<q( t) q*( t+r )>  =  <S>p2(x) 
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(2-21) 
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The  correlation  function  for  the  bistatic  case  is 

<q{t)q*(t+T)>  =  <  p1(t)p2(t)pi*(t+T)p2*(t+T)> 

=  <  Pi(t)pi*(t+r)><p2(t)p2*(t+T)>  {2-22) 

when  it  is  assumed  that  pi(t)  and  P2(t)  are  independent.  In  this  case, 

<q( t )q*( t +t  )>  =  <S>pj ( t  )  p2 ( t)  (2-23) 

In  qeneral  for  a  Gaussian  PSD, 

<q(t )q*(t+T)>  =  <S>  exp[-(T/TQ)2]  (2-24) 

where  i0  is  the  decorrelation  time  of  the  received  voltage  and 

one-way 

Ti //2  monostatic 

(2-25) 

bistatic  (*1,1*11,2) 

Ti. i+Ti.i 

Ti,i/^  bistatic  (t^j  =  t1<2) 

For  bistatic  propagation,  Tj^i  =  1  or  2!,  is  the  one-way  decorrelation 
time  of  the  ith  path.  It  can  be  seen  for! the  bistatic  case,  that  if  one 
of  the  paths  has  a  small  decorrelation  time  relative  to  the  ether  one-way 
path,  the  resulting  decorrelation  time  of  ^he  received  signal  corresponds 
to  the  smaller  decorrelation  time.  If  both  one  way  paths  have  the  same 
decorrelat ion  time,  the  resulting  t0  corresponds  to  the  ironostatic  case 
with  a  decorrelation  time  of  \\H2.  In  uorh  two-way  cases,  it  is 
apparent  that  the  decorrelation  time  of  the  received  signal  is  smaller 
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than  the  decorrelation  time  of  the  one-way  propagation  path.  For 
arbitrary  power  spectral  densities  and  propagation  geometries,  the 
decorrelation  time  of  the  received  voltage  is  related  to  the  one-way 
decorrelation  time  as 


T 

0 


I 

I 


T1 

6x^/2 

/  7"2 

6t1, 1 T 1 , 2/' T 1 ,  1 


+ 


2 

>1,2 


one  way 
monost  at  ic 

bistat  ic 


(2-26) 


The  scale  factor  5  differs  from  unity  in  general  if  the  PSD  is  r.ot 
Gaussian. 


2.4  FIRST  ORDER  STATISTICS  FOR  DUAL  CHANNEL  COMMUNICATION  SYSTEMS 

Dual  channel  communication  systems  may  combine  the  amplitudes  of 
two  independent  channels  in  order  to  mitigate  the  effects  of  propagation 
fading.  Two  possible  algorithms  for  combining  the  amplitudes  of  the  two 
channels  into  the  amplitude  z  used  for  demodulation  are 

(  rl  +  r2  summation  algorithm 
2=[  (2-27) 

|max(ri,r2)  maximum  algorithm 

where  rj,  i=l  or  2,  is  the  amplitude  of  the  received  voltage  on  the  i  th 
channel  and  the  function  max(x,y)  is  equal  to  the  large  of  x  and  y.  If  it 
is  assumed  that  the  two  channels  utilize  propagation  paths  of  sufficient 
spatial  separation  that  the  propagation  effects  on  the  received  voltages 
are  independent.  Then  in  the  strong  scattering  limit,  r!  and  r2  are 
independent,  Rayleigh  distributed  random  variables.  It  will  also  be 

assumed  that  the  transmitted  power  is  equal  in  the  two  channels  so  that 

2  2 

<rj>  =  <r2>. 
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for  the  summation  algorithm,  the  cumulative  distribution  of  z 
has  been  given  by  Marcum5  from  which  the  cumulative  distribution  of  the 
combined  power  (S  -  z2)  is  obtained  as 

F(S)  =  l-exp[ -S/2o2]  -  A  S/4 a2  e'S/4°2  erf (/s/4o2 )  (2-28) 

where  2a2  is  the  mean  power  in  each  of  the  channels  and  erf  is  the  error 
function.  The  mean  power  of  the  combined  amplitude  is 

<S>  =  <z2>  =  <r\>  +  2<r1Xr2>  +  <r2>  =  (4  +  n)a2  .  (2-29) 

The  cumulative  distribution  of  z  for  the  maximum  algorithm  is 

F(z)  *  (1  -  expf-z2 /2o2 )]2  (2-30) 

where  the  term  inside  the  brackets  is  the  probability  that  one  of  the 
Rayleigh  distributed  amplitudes  is  less  than  or  equal  to  z  and  the  square 
gives  the  probability  that  both  rx  and  r2  are  less  than  or  equal  to  z. 

The  mean  received  power  is  then  easily  calculated  as 

<S>  =  <z2>  =  /  z2  dF(z)  =  3o2  '  (2-31) 

o 

and  the  cumulative  distribution  of  the  combined  power  is 


F(S)  =  (1  -  exp[ -3S/2< S>] } 2  . 


(2-32) 


The  probability  density  function  of  the  combined  power  is  obtained  by 
differentiating  Equation  2-32  qiving 

f(S)  *  (3/<S>)  exp(-3S/2<S>)  [l  -  exp( -3S/2<S>]  (2-33) 

The  median  combined  power  for  the  d.:3l  channel  communication 

systems  is 


<4u0/(4+if)  =  -0.81dB  summation  algorithm 

M/<S>  =  J  '  (2-34) 

(2  tn[/2/(/2-l)  ]  *  -0.87dB  maximum  algorithm 

where  u0  is  the  solution  of  the  equation 

e  0  +  /iruQ  e  0  erf(/uQ)  =  1/2  (2-35) 

which  yields 

uQ  *  1.482227  .  (2-36) 

When  compared  with  the  single  channel  median  power  of  1.6  dB  below  the 
mean  power,  these  results  indicate  that  the  dual  channel  system  does  in 
fact  mitigate  the  effects  of  scintillation. 

The  cumulative  distributions  for  the  combined  power  are  plotted 
on  probability  paper  in  Fiqure  2  for  the  two  dual  channel  algorithms  along 
with  the  cumulative  distr ibut w  oT  the  received  power  of  the  one-way 
communication  link  or  single  channel  case.  The  near  coincidence  of  the 
dual  channel  curves  is  a  consequence  of  the  fact  that  because  the  ampli¬ 
tudes  from  the  two  channels  are  independent  and  Rayleigh  distributed,  it 
is  probable  that  one  amplitude  is  much  smaller  than  the  other  and  hence 
that  the  larger  amplitude  is  nearly  equal  to  the  sum  of  the  two  ampli¬ 
tudes.  It  is  apparent  from  the  figure  that  the  dual  channel  system  has  a 
lower  probability  of  deep  fades  than  does  the  single  channel  system.  For 
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Figure  2.  Cumulative  distribution  of  the  signal  power  used  for  demodula¬ 
tion  in  single  and  dual  channel  communication  links. 

example,  the  dual  channel  system  has  a  probability  of  about  3  x  10*1*  of 
fades  of  20dB  or  greater  while  this  probability  is  10*2  for  the  single 
channel  system. 

2.5  AMPLITUDE  MOMENTS 

The  expected  values  of  the  first  4  moments  of  the  received 
amplitude  z,  where 

2  "  |q(t)|  ,  -  (2-36) 
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for  one-way,  monostatic,  a;.d  bistatic  cases,  and  for  the  combined  ampli¬ 
tude  defined  in  Equation  2-26  for  the  dual  channel  communications  systems 
are  needed  for  comparison  with  the  simulation  results.  These  are  given  in 
Table  1  and  can  be  calculated  directly  from  the  moments  of  a  Rayleigh 
distribution  for  all  but  dual  channel  maximum  algorithm.  In  the  latter 
case.  Equation  2-29  is  used  to  calculate  the  moments  of  z. 


Table  I.  Moments  of  the  amplitude  of  the  received  voltage. 
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SECTION  3 

BINARY  ERROR  PROBABILITIES  FOR  TWO-WAY 
COMMUNICATION  LINKS  UNDER  SLOW  RAYLEIGH 


AND  DUAL  CHANNEL 
FADING  CONDITIONS 


Rinary  error  probabilities  of  two-way  and  dual  channel  communi¬ 
cation  links  are  derived  in  this  section  for  noncoherent  M-ary  frequency 
shift  keyinq  (M-ary  FSK).  coherent  phase-shift  keying  (PSK),  differen¬ 
tially  encoded  coherent  PSK  (APSK),  and  different ial  ly  coherent  binary  PSK 
(0BPSK)  under  slow  Rayleigh  fadinq  conditions.  Perfect  phase  tracking 
will  be  assumed  for  the  PSK  demodulation  techniques,  *n  the  two-way 
qeometry,  a  siqnal  is  transmitted  throuqh  an  ionospheric  channel  to  a 
transponder  which  then  retransmits  the  siqnal  throuqh  an  independent 
ionospheric  channel  to  a  receiver  which  is  not  colocated  with  the  original 
transmitter.  It  is  assumed  that  the  transponder  is  an  ideal  linear  ampli¬ 
fier  that  is  noiseless  and  that  outputs  a  siqnal  with  identical  statistics 
to  the  input  siqnal.  For  a  dual  channel  communication  link,  independent 
siqnals  are  obtained  by  the  use  of  spatial  diversity  (e.g.,  separated 
receiver  antennas  or  transmitter  locations)  or  by  other  diversity  tech¬ 
niques  such  as  multiple  transmission  frequencies.  The  results  of  the 
previous  section  showed  that  the  distributions  of  the  combined  signals  for 
the  two  algorithms  considered  are  essentially  identical.  Therefore,  the 
binary  error  prohabi  1  it ies  will  be  calculated  only  for  the  maximum  algo¬ 
rithm  which  has  the  simpler  mathematical  form  for  the  distribution  of 
combined  power. 


The  probability  of  bit  error  in  an  additive  white  Gaussian  noise 
(AWGN)  channel  for  a  constant  energy  siqnal  and  for  M-ary  FSK  modulation 


P  (E  /N  )  =  1  f 

be  b  o  2(M-1 ) 


M 

I 

k=2 


■1) 


(") 


expl-(^J.)(loq  M)lH] 
k  2  N0 


(3-1) 


25 


Upon  changing  variables  (u  =  2/Eb/<Eb>) ,  the  integral  reduces  to 

I  =  /  u  exp[  -a<Eb/NQ>u2  /4]  Kq(u)  du  .  (3-1 

This  integral  is  given  by  Gradshteyn  and  Ryzhik7  in  terms  of  one  of 
Whittaker's  function  which  can  be  further  reduced  using  the  mathematical 
properties  of  confluent  hypergeometric  functions  from  Abramowitz  and 
Stegun8  to  give 


exp[l/a<E./N  >]  E1[l/a<E./N  >] 


a<E. /N  > 
b  o 


(3-7) 


where  Ei  is  the  exponential  integral. 


The  average  probability  of  bit  error  is  then 

<P  >  = _ I _  y 

be  2(M-1)( log2M)<E  /N  >  1  '  '  \k /  k-1 


b'  o  k =2 
k/( k-1 ) 

( 1 og2  M)  <Eb/NQ 


(3-8) 


>]  E^[ 


k/(k-l) 

(log2M)  <Eb/No> 


for  M-ary  FSK  demodulation  and  is 


s  _  exp[  l/<Eb/Mo>]  E,[l/<Eb/No>] 

y  *  ~  1  r  ^  jm  ,»i|"  ■— 


(3-9) 


for  DBDSK  demodulation.  The  mean  probability  of  bit  error  cannot  to  the 
authors  knowledge  be  obtained  in  closed  form  for  ideal  coherent  PSK  or 
aPSK  modulation.  However,  by  using  the  same  change  of  variables  that  was 
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used  to  derive  Equation  3-6,  the  mean  probability  of  bit  error  for  these 
cases  can  be  rewritten  as 


<P.  >  = 

be 


(1/2)  /  u  Kq(u)  erfc[/Eb/NQ>  u/2]  du 


ideal  PSK 


(3-10) 


/  u  K  (u)  erfc[/<E./N  >  u/2]  {1  -  (1/2)  erfc[/<E  /N  >  u/2]}  du 


ideal  aPSK 

Both  integrals  of  Equation  3-10  nave  well  behaved  integrands  and  are 
easily  evaluated  numerically. 

For  dual  channel  links  under  slow  Rayleigh  fading  conditions, 
'.he  average  bit  error  probability  with  M-ary  FSK  demodulation  can  be 
evaluated  directly  to  give 


<P  >  .  - L_  l  (.!)«  (M) 

be  2(M-1)  Vk/ 


(3-11) 


[3+(k-l)(log2M)<Eb/No>/k]  [ 3+2(k-l) ( log2M)<Eb/No>/k] 


Similarly,  with  DBPSK  demodulation. 


<p  >  =  _ lit _ 

be  (3+<Eb/No>)(3+2<Eb/No>)  * 


(3-12) 


The  evaluation  of  Equation  3-4  for  ideal  PSK  and  ideal  aPSK 
involves  integrals  of  the  form 


I  =  /  erfc(/x)  e'aX  dx  =  (2/a)  ( 1-1//TCT) 
0 


(3-13) 


•*  # 

4 


and 


I  *  /  erfc2 (/x)  e“aX  dx 
o 


(3-14) 


/1+a 


where  the  first  integral  is  given  by  Gradshteyn  and  Ryzhik  and  the  second 
integral  can  be  integrated  by  parts  to  reduce  it  to  a  form  given  by 
Gradshteyn  and  Ryzhik.  The  mean  probability  of  bit  error  for  ideal  PSK  | 

i 

then  becomes 


<V  ■  i - & 


1 


/l+3/<Efc/N0>  /l+3/(2<Eb/No>) 


and  for  ideal  a PSK  it  becomes 


<P 


j  [ l-( 2/w )  tan-1  /l+3/<Eb/N0>] 
>  =  _  +  "  ■  — ■  •  ■ 


be  "  7 


/l+3/<E./N  > 
b  o 


(3-16) 


2[l-(2/n)  tan-1  /l+3/2<E. /N  >] 

DO 


/l+3/2<E./N  > 
b  o 


The  probabilities  of  bit  error  for  one-way  or  single  channel 
communication  links,  for  two-way  links,  and  for  dual  channel  links  under 
slow  Rayleigh  fading  conditions  are  shown  in  Figures  3  through  8  for  non¬ 
coherent  M-ary  FSK  with  M  equal  to  2,  4,  and  &,  ideal  coherent  PSK,  a^SK, 
and  D8PSK  demodulation  respectively.  The  probability  of  bit  error  for  an 
AWGN  channel  is  also  shown  for  each  of  the  techniques.  Note  that  an  ideal 
coherent  PSK  dual  channel  system  is  not  achievable  in  practice  because  of 
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Figure  3.  Theoretical  performance  of  noncoherent  binary  FSK  demodulation 
in  an  AMGN  channel  and  under  slow  Rayleigh  fading  conditions. 


Figure  4.  Theoretical  performance  of  noncoherent  quaternary  FSK  demodula¬ 
tion  in  an  AWGN  channel  and  under  slow  Rayleigh  fading 
conditions. 
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Figure  5-  Theoretical  performance  of  noncoherent  8-ary  FSK  demodulation 
in  an  AWGN  channel  and  under  slow  Rayleigh  fading  conditions. 


Figure  6.  Theoretical  performance  of  ideal  coherent  PSK  demodulation  in 
an  AWGN  channel  and  under  slow  Rayleigh  fading  conditions. 
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Figure  7.  Theoretical  performance  of  ideal  coherent  &PSK  demodulation 
in  an  AWGN  channel  and  under  slow  Rayleigh  fading  conditions. 


Figure  8.  Theoretical  performance  of  differentially  coherent  binary  PSK 
demodulation  in  an  AWGN  channel  and  under  slow  Rayleigh  fading 
conditions. 
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the  k  phase  ambiguity.  As  can  be  seen  from  the  figures,  the  mean  proba¬ 
bility  of  bit  error  for  a  given  bit  energy-to-noise  density  ratio  is 
higher  for  a  two-way  link  than  it  is  for  a  one-way  link  because  of  the 
increased  probability  of  occurrence  of  deep  fades  in  the  two-way  geom¬ 
etry.  However,  the  mean  probability  of  bit  error  for  a  given  bit  energy- 
to-noise  density  ratio  is  lower  for  a  dual  channel  link  than  it  is  for  a 
single  channel  link  because  of  the  decreased  probability  of  occurrence  of 
deep  fades  in  the  dual  channel  link.  At  a  bit  error  probability  of  10“3, 
the  two-way  link  requires  between  8.5  dB  and  9.5  dB  more  bit  enerqy-to- 
noise  density  than  does  the  one-way  link  for  all  demodulation  techniques 
considered.  Conversely,  the  dual  channel  link  requires  between  9.5  dB  and 
11  dB  less  bit  energy-to-noise  density  than  does  the  single  channel  link 
to  achieve  a  10-3  bit  error  probability. 


SECTION  4 

LEVEL  CROSSING  PROBLEM 


The  main  point  of  this  study  is  a  statistical  characterization 
of  the  2ero  crossinqs  of  a  random  process  z'(t)  defined  as 

z'(t)  =  z(t)  -  *  (4-1) 

where  i  is  an  arbitrary  amplitude  level.  This  problem  was  first  treated 
in  detail  by  Rice9  for  the  case  that  z(t)  is  a  normally  distributed  random 
variable.  The  purpose  of  this  section  is  to  extend  the  results  of  Rice  to 
the  cases  where  z(t)  is  the  amplitude  of  q(t)  defined  in  Equation  2-9. 

For  one-way,  monostatic,  or  bistatic  cases,  an  interesting  statistic  of 
the  zero  crossing  problem  that  can  be  calculated  in  closed  form  is  the 
mean  number  <N( A ,T)>  of  crossings  of  the  level  *  in  the  time  interval  T. 
From  the  mean  number  of  level  crossings,  the  mean  duration  and  separation 
of  fades  and  f’ares  can  also  be  calculated  in  closed  fo>m. 

It  can  be  shown  that9 
t  '+T  • 

<N(*,T)>  =  /  dt  /  dz  | z |  f(*,z;t)  (4-2) 

{  • 

where  f(z,z;t)  is  the  joint  probability  density  function  of  z  and  its  time 
derivative  z  at  time  t.  For  the  stationary  processes  considered  here, 
f(z,z;t)  is  independent  of  t  and 

OB 

•  <N(*,T)>  =  T  /  |z|f(*,z)dz  .  (4-3) 

»O0 
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4.1 


ONE-WAY  PROPAGATION 


The  first  step  in  evaluating  Equation  4-3  is  to  calculate  the  joint  proba¬ 
bility  density  function  of  the  components  of  the  one-way  received  voltage 
and  their  time  derivatives.  The  joint  probability  density  function  of  the 
in-phase  and  the  quadrature  phase  components  of  the  received  one-way  volt¬ 
age,  denoted  by  x  and  y  respectively,  is  given  by  Equation  2-2.  Rice  has 

shown  that  if  x  and  y  are  normally  distributed  and  stationary  with  zero 

•  • 

means,  then  x  and  y  are  normally  distributed  and  stationary  with  zero 
means  and  variances  given  by 


o?  =  a?  =  -  a2p(o)  =  2o2A2/t1  (4-4) 

x  y 


where  p(t)  is  the  seconH  time  derivative  of  the  autocorrelation  function 
of  the  one  way  voltaqe.  In  addition,  x  and  y  are  independent  of  both  x 
and  y.  For  arbitrary  power  spectral  densities,  it  is  then  required  that 
p(o)  be  finite.  The  implications  of  this  requirement  are  discussed  in 
more  detail  in  Appendix  A.  The  scalinq  factor  A2  is  unity  for  a  Gaussian 
PSD  and  is  calculated  for  a  f-M  PSD  in  the  Appendix.  Thus  the  joint 
probability  density  function  of  x,  x,  y,  and  y  is 


f(x,x,y,y) 


exp{-[x2+y2+  Ti(x2+y2)/?A2]/2a?} 
(  2jt  o2 ) 2  (2A2/t2)  “ 


(4-5) 


The  joint  probability  density  function  of  the  amplitude  r  of  the 
one-way  received  voltaae  and  of  r  can  now  be  obtained  from  Equation  4-5  by 
chanqinq  variables.  Lettinq 

x  =  r  cose  (4-6) 
y  *  r  sine 
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where  0  is  the  phase  of  the  received  voltage,  the  time  derivatives  of  x 
and  y  become 


so  that 


x  *  r  cos9  -  r9  sin6 
y  *  r  sin6  +  r0  cos8 


x2  +  y2  =  r2 


(4-7) 


i 

<  ■.  'S- 

\  •  V 


x2  +  y'-  =  r 2  +  r202 


(4-8) 


•  m 

Equation  4-5  can  be  then  be  rewritten  jin  terms  of  r,  r,  0,  and  0  as 

f(x,x,y,y)  dxdrdydy  *  f(r,r,0,*0)  jdet ( J)  | drdrdOdO  (4-9) 

where  the  determinant  of  the  Jacobian  J  of  the  coordinate  transformation 
is  ! 


det(J)  =  det 


3x/3r 

3x/3ij- 

3y/3r 

3y/3r" 

3x/3r 

3x/3f 

3y/3r 

3y/3r 

3x/30 

3x/30 

3y/39 

3y/30 

3x/30 

3x/30 

3y/39 

3y/30j 

cos8 

-esin9 

Sin9 

0COS0 

0 

cos9 

0 

sin0 

-rsin8 

-rsin9- 

-r9cos9 

rcos9 

fcos0 

0 

-rsin9 

0 

rcos9 

(4-10) 


w%  ;x 

i  -A  ■ 
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j  V-  i 

I  \  .y ‘  t 
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i  •  ^ 
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•  V  ’  \ 
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The  joint  probability  density  function  of  r  and  r  is  qiven  by  inteqratinq 
Equation  4-9  over  0  and  0: 

2  IT  «> 

f(r,r)drdr  =  J  (10  J  d9  f(r,r,0,0)  drd r  (4-11) 

o  -® 

r2  pxp[-(r2+t2  r2/2A2)/2cr2  jdrdr 
(2*o2)2  (2A2/t2) 

*  /  d9  /  d0  exp(-t2r 2&2/4A2o2) 
o  -«•  1 

r  exp(-r2 /2a2)  exp(-T2r2/4A2o2) 

=  _  dr  •  - - L... — _  dr 

/4itA2o2/t2 

1 

It  can  bp  seen  from  Fquation  4-11  that  the  probability  density  function  of 
r  is  Raylpiqh;  thp  probability  density  function  of  r  is  Gaussian  with  zero 

moan  and  2A2o2/t1  standard  deviation;  and  because  f(r,r)  is  separable  into 

•  • 

a  function  of  r  times  a  function  of  r,  r  and  r  are  independent. 

Thp  mean  number  of  crossinq  in  the  time  interval  T  for  the 
onp-way  rasp  can  now  be  easily  evaluated  from  Equations  4-3  and  4-11. 
dsinq  Equation  2-1?  to  write  o  in  terms  of  the  mean  received  power,  the 
mean  number  of  crossinq  becomes 

<N(L,T)>  =  a(T/to)  e'L/<S>  (4-12) 

where  l  =  t2  is  the  power  of  the  level  l  and  where  t0  is  equal  to  tj  in 
this  case.  The  effect  of  arbitrary  power  spectral  densities  is  then  to 
scale  thp  mean  number  of  IpvpI  crossinqs  in  the  one-way  propaqation  by  the 
quantity  A.  Notinq  that  two  crossinqs  are  required  for  a  fade  (flare) 
below  (above)  the  level  l,  the  expected  number  of  fades  (flares)  per  unit 
time  is  >N(L,T)>/2T. 
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4.2 


MONOSTATIC  PROPAGATION 


For  the  monostatic  case,  the  amplitude  z  of  the  received  voltage 
is  equal  to  r2  so  that  x  and  y  as  functions  of  z  and  the  phase  0  are 


1  /2 

x  =  z  cose 

1  /2  . 

y  ■  z  sine 


(4-13) 


Using  this  transformation,  the  calculation  of  f(z,z)  is  then  identical  to 
the  one-way  case.  In  this  case,  the  time  derivatives  or  x  and  y  become 


so  that 


•  “1  /?  •  ,  1/2  • 

x  =  z  z  cose/2  -  z  e  sine 

•  •l/2«  ,  1/2  • 

y  =  2  z  sine/2  +  z  0  cose 


x2+y2  *  z 

x2  +V2  =  z2/4z  ♦  ze2  . 


(4-14) 


(4-15) 


Again  the  change  of  variables  requires  the  evaluation  cf  the  determinant 
of  the  Jacobian,  defined  in  Equation  4-10  by  replacing  r  with  z,  which 
gives 


l'i/}  cost/?  l  cost/t  -  l *,/l  t  tint/?  r,/f  tint/?  i  tint/*  ♦  j‘,/?  t  cott/? 


dft(J)  •  dot  I 


cost/? 


f',/*  s  ln»/r 


-I1'*  slot  l"1'1  i  tint/?  -  i  com  i1'1  cost  i'”1  i  cost/?  -  tl,t  i  tint 


•tlf*  stno 


l'h  cost 


(4-16) 


**  r 

•  ■  1 
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The  joint  probability  density  function  of  z  and  z  can  now  be  written  as 


f(z,z)dzdz  *  (1/4)  /  de  /  de  f(z,z,e ,e )dzdz 
o  -» 


exp[-(z+t2  z2 /8a2  z) /2o2  ]  ^  "  .  2  2  , 

- —  £  ®i  d8  exp(-’>  /4‘ 0  )d2d2 


(4-17) 


=  (t^ /8Ao3/iTz)  exp[-(  z  +  T2*z2/8A2z)/2o2]dzdz  . 


Because  Equation  (4-17)  is  not  separable  into  a  function  of  z  times  a 
function  of  z,  z  and  z  are  statistically  dependent. 

The  mean  number  of  crossings  for  the  monostatic  case  of  the 
level  L  per  time  interval  T  is  easily  evaluated  from  Equations  4-3  and 
4-17  as 


<N( L ,T)>  =  «a(T/tq)  [32L/h2<S>]  exp[ -/2L/<S>] 


(4-18) 


where  t  is  equal  to  5t.//2  for  two-way  propagation  and  a  is  again  written 
o 

in  terms  of  the  mean  power  using  Equation  2-12. 

4.3  BISTATIC  PROPAGATION 

For  the  bistatic  propagation  case,  the  amplitude  z  of  the 
received  voltage  q(t)  from  Equation  2-9  can  be  written  as 


where  s  and  t  are  the  amplitudes  of  the  two  independent  bistatic  path 
voltages  px(t)  and  p2(t).  The  joint  probability  density  function  of  s,  s, 
t,  and  t  is  then 

f(s,s,t,t;  -  ■  vi, a)  ( 4-20) 


•  • 

where  f(s,s)  and  f(t,t)  are  both  given  by  Equation  4-11  assuming  that  each 
path  has  the  same  PSD  and  decorrelation  time.  The  more  general  problem 
where  the  PSD  and/or  the  decorrelation  time  of  the  two  paths  differ  is 
considered  in  Appendix  B. 


In  this  case  it  is  convenient  to  use  Equation  4-19  to  write 
•  •  •  • 
f(s,s,t,t)  in  terms  of  a  new  set  of  variables  z,  z,  s,  and  s.  The 

Jacobian  of  this  transformation  is 


3S/3Z 

3S/3Z 

3  t/3  Z 

3  t/3  Z 

J  = 

3  s/3  Z 

3  s/3  Z 

• 

3t/3Z 

3 1/3  Z 

• 

(4-21) 

3  s/3  S 

3  S/3  S 

3  t/3  s 

3 1/3  S 

• 

3  s/3  S 

•  • 

3  s/3  S 

3  t/3  S 

•  • 

3  t/3  S 

where  the  old  variables 

• 

(s,  s, 

• 

t,  and  t) 

in  terms  of 

the  new  variables 

•  • 

(s,  s,  z,  and  z)  are 


s  =  s 

•  • 

s  *  s 

t  =  z/s 

t  -  i/s  -  z;/s2 

so  that  the  determinant  of  the  Jacobian  becomes 
det(J)  =  s-2 


(4-22) 


(4-23) 


« 


I 
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The  joint  probability  density  function  of  z  and  z  is  then  obtained  by 

•  •  • 

integrating  f(s,s,z,z)  over  s  and  s  as 


f(z,z)  =  (t2/4jta2o6)  /  ds  /  ds  zs-2 
1  o  — 

(4-24) 

x  exp{-[s2  +  z2/s2  +  t2(s2  +  *z2s-2  +  z2s2s“4  -  2zzs/s3)/2Az]/2a2} 

The  s  integral  can  be  performed  in  closed  form  to  give 

|  f(z,z)  =  (tj z/2/TAa5 )  /  ds  (s^+z2)"1^2 

o 

x  exp(-[s2  +  zVs2  +  t2  s2  z2  ( s4  +z2 ) -1  /2a2  ]  /2a2 }  (4-25) 

|  ' 

The  mean  number  of  crossings  of  the  level  L  requires  an  integral 

Of  Equation  4-25  over  z  which  can  be  obtained  in  closed  form.  Upon  chang¬ 
ing  variables  in  the  remaining  integral  to  u  =  s/a;  using  Equation  2-12  to 
write  o  in  terms  of  the  mean  power;  and  noting  that  t0  is  equal  to  6t1//2 
in  tne  bistatic  geometry  when  the  two  paths  have  the  same  decorrelation 
t|ime;  the  mean  number  of  crossings  becomes 

•  | 

<N(L,T):  *  6A(T/T0)/8L/<S>ir  /  /l+4L/<S>u4  exp[ -(l+4L/<S>u4  )u2/2]  du 

o 

(4-26) 

The  remaining  integral  is  easy  to  evaluate  using  numerical  techniques 
because  the  exponential  term  in  the  integrand  is  exp(-2L/<S>u2 )  for  u  «  1 
and  is  exp(-u2/2)  for  u  »  1.  Hence  the  integrand  goes  exponentially  to 
zero  at  both  limits  of  the  integral. 


41 


4.4 


MEAN  DURATION  AND  SEPARATION  OF  FADES  AND  FLARES 


The  mean  numbers  of  level  crossings  per  received  decorrelation 
time  are  plotted  in  Figure  9  versus  the  ratio  L/<S>  in  dB  for  the  one-way, 
monostatic,  and  bistatic  propagation  cases  and  for  a  Gaussian  PSD  for 
which  a  and  a  are  both  unity.  These  results  show  the  higher  frequency  of 
deep  fades  in  the  two-way  monostatic  or  bistatic  propagation  cases  versus 
the  one-way  propagation  case.  They  also  show  that  for  levels  between  +5 
dB  and  about  -12  dB  relative  to  the  mean  power,  the  bistatic  case  has  a 
higher  frequency  of  level  crossings  than  does  the  monostatic  case. 


LE  VEL/<PO  WER  >  (dB) 


Figure  9.  Mean  number  of  level  crossings  in  one  decorrelation  time  t0. 


The  mean  duration  <tg>  of  fades  or  flares  can  be  obtained 
directly  from  the  mean  number  of  level  crossings  by  noting  that  10 


(1/2)  <*U,T0»  <t0ifade>/.o  -  F(L) 
(1/2)  <»(L„0)>  <tD>flare>/T0  ■  1  -  F(L) 


(4-2?) 


where  F(L)  is  the  probability  that  the  received  power  is  less  than  or 
equal  to  L  (Equation  2-11)  and  (1/2)  <N(L,t  )>/t  is  the  mean  number  of 

O  9 

fades  or  flares  per  unit  time.  In  addition,  the  mean  separation  <t$> 
fades  or  flares  is 


<tS,fade>  ’  <tS, flare  >  *  V<N<L>To»  «-28> 

or  just  twice  the  mean  interval  between  level  crossings. 

The  scaling  of  the  mean  duration  and  separation  of  fades  or 
flares  from  the  Gaussian  values  <tg>g  and  <t$>g  is 


D  "  6A  D  G 


<tS>  =  6A 


(4-29) 


where  again  6  is  unity  for  one-way  propagation.  Under  Rayleigh  fading 
conditions,  the  PSD  will  be  bounded  by  the  Gaussian  form  in  the  strong 
scatter,  fast  fading  limit  with  1/6A  =  1  and  by  an  f-3  PSD  in  the  slow 
fading  limit.  An  f-4  spectrum  lies  between  these  two  limits  and  is  of 
particular  interest  since  it  is  frequently  used  for  system  specification 
and  testing.11  For  an  f-4  PSD,  the  values  of  1/5A  that  are  calculated  in 
Appendix  A  are  0.66  for  one-way  propagation  and  0.74  for  two-way  propaga¬ 
tion  when  both  of  the  one-way  paths  have  the  same  decorrelation  time. 

The  mean  duration  of  fades  below  an  arbitrary  level  L  for  a 
Gaussian  PSD  is  shown  in  Figure  10  and  the  mean  duration  of  flares  above 
the  level  L  for  a  Gaussian  PSD  is  shown  in  Figure  11.  The  monostatic  case 
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RAYLEIGH  FADING 
GAUSSIAN  PSD 


LEVEL/ <P0  WER>  (dB) 

Figure  10.  Mean  duration  cf  fades  below  an  arbitrary  level 


RAYLEICH  FADING 

GAUSSIAN  PSD 

LEVEL/ <P0  WER>  (dB) 

Figure  11.  Mean  duration  of  flares  above  an  arbitrary  level 


has  the  longest  mean  duration  of  fades  for  any  level  L  and  the  one-way 
case  has  the  shortest  mean  duration  of  fades.  However,  above  levels  of 
about  -3  dB  relative  to  the  mean  power,  the  monostatic  case  also  has  the 
longest  mean  duration  of  flares  while  the  one-way  case  has  the  shortest 
mean  duration  of  flares.  The  situation  is  reversed  at  levels  below  -3 
dB.  That  is,  the  one-way  case  has  the  longest  mean  time  above  levels 
which  are  3  dB  or  more  below  the  mean  power.  The  bistatic  mean  duration 
of  fades  or  flares  lies  between  the  other  two  cases  except  in  the  cross¬ 
over  region  (-10  dB  to  -3  dB)  of  the  mean  ouration  of  flares. 

The  mean  separation  of  fades  or  flares  for  the  Gaussian  PSD  is 
shown  in  Figure  12.  The  monostatic  case  has  shortest  mean  separation  of 
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fades  or  flares  for  levels  11  dB  or  more  below  the  mean  power  and  for 
levels  8  dB  or  more  above  the  mean  power.  For  levels  between  -9  dB  and  2 
dB  relative  to  the  mean  power,  the  one-way  case  has  the  shortest  mean 
separation  of  fades  or  flares.  Except  where  the  curves  cross  over,  the 
bistatic  mean  separation  of  fades  or  flares  lies  between  the  other  two 
cases. 


Using  these  results,  qualitative  information  can  be  obtained  on 
the  difference  between  mean  performance  effects  of  a  Rayleigh  fading 
environment  on  one-way  versus  two-way  communication  links  or  monostatic 
versus  bistatic  space  based  radars.  In  the  communication  link  case,  a 
comparison  of  the  one-way  and  bistatic  propagation  cases  shows  that  the 
mean  duration  of  fades  is  longer  for  the  two-way  link;  that  the  mean  dura¬ 
tion  of  flares  above  the  mean  power  is  longer  for  the  two-way  link;  and 
that  for  levels  between  -9  dB  and  2  dB  relative  to  the  mean  power,  the 
one-way  link  has  shorter  mean  separations  between  fades  or  flares  while 
the  one-way  link  has  longer  mean  separations  between  fades  or  flares  for 
levels  outside  of  the  -9  dB  to  2  dB  range.  By  comparing  the  monostatic 
and  bistatic  propagation  cases  it  can  be  seen  that  the  mean  duration  of 
fades  is  longer  for  the  monostatic  radar;  that  the  mean  duration  of  flares 
above  the  mean  power  is  also  longer  for  the  monostatic  radar;  and  that  the 
bistatic  radar  has  shorter  mean  separations  of  fades  or  flares  for  levels 
between  -11  dB  and  3  dB  relative  to  the  mean  power  while  the  monostatic 
radar  has  longer  mean  separations  of  fades  or  flares  outside  of  the  -11  dB 
to  3  dB  range. 
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SECTION  5 

GENERATION  OF  THE  RECEIVED  VOLTAGES  AND  DATA  GATHERING 


The  probability  density  functions  of  the  durations  and  separa¬ 
tions  of  fades  or  flares  are  unobtainable  by  analytic  methods  for  even  the 
simplest  case  of  a  one-way  communicat ion  link.  However,  the  desired 
statistics  are  easily  calculated  using  Monte-Carlo  methods  by  generating 
either  one  or  two  random  sequences  of  the  one-way  voltage  p(t)  with  Gaus¬ 
sian  correlation  and  then  using  Equation  2-9  to  generate  the  sequences  of 
the  received  voltage  q(t).  It  is  then  an  easy  matter  to  find  the  statis¬ 
tics  of  the  duration  separation  of  fades  or  flares  relative  to  a  given 
level  that  exist  in  the  amplitudes  z(t).  The  sequences  generated  are 
limited  to  be  less  than  or  equal  to  8192  samples  because  of  the  fast 
Fourier  transform  routine  that  is  utilized.  Thus  this  process  is  repeated 
multiple  times  in  order  to  collect  sufficient  data. 

The  method  of  generating  the  sequences  is  described  in  detail 
elsewhere12  and  is  briefly  reviewed  here.  The  desired  Gaussian  autocor¬ 
relation  is  equivalent  to  a  Gaussian  PSD  S(f)  where 


S(f)  =  /  p(t)  exp(-i2*ft)dt 

=  r  /ir  exp(-ir2T2f2  ). 

1  1 


(5-1) 


The  voltage  samples  are  first  generated  in  frequency  space  using  the 
discrete  equivalent  Sn  of  the  PSD  S(f)  and  are  then  discrete  Fourier 


t  * 

*»  I 


47 


transformed  to  the  time  domain.  A  total  of  N=2m  (m  is  a  positive 
integer  less  than  a  equal  to  13)  voltage  samples  are  generated  covering  a 
total  time  of 

T  =  NAt  .  (5-2) 

If  there  are  n  samples  per  one-way  decorrelation  time,  then 

At.  =  ij/n  .  (5-3) 

Because  the  spectral  samples  of  a  stationary  random  process  are  independ¬ 
ent,  N  complex  voltages  cj  (j=0,l,...,  N-l)  are  generated  with  Rayleigh 
distributed  amplitudes  and  unifrmly  distributed  phases.  Thus 

5  j  .  /-2a24n(u1  j)  exp(  i2n u2  j) ,  j  *  0,1 . N-l  (5-4) 

where  Uj  j  and  Uj j  are  independent,  uniformly  distributed  random  varia¬ 
bles  on  the  interval  [0,1).  The  real  and  imaginary  parts  of  are 
independent,  normally  distributed  random  variables  with  zero  means  and 
o2  variances.  The  parameter  o  is  set  by  requiring  that  the  average 
received  power  be  unity  in  all  cases.  The  frequency  samples  are  then 
weighted  symmetrically  by  the  discrete  voltage  spectral  density  /sT  where 

J 

Sj  ■  Af  S[(j-N/2)Af],  j  =  0,1,. ...N-l,  (5-5) 

and  where  the  frequency  sample  spacing  is 

Af  =  1/T  *  1/NAt  =  n/Nti  .  (5-6) 
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Hence,  the  frequency  voltaqes  xj  are  given  by 

Xj  =  Cj  (n/w  /N) 1 exp[-*2n2(j-N/2)2/(2N2)] 


j  =  0,1 .... ,N-1 


The  one-way  received  voltage  p.  is  the  discrete  Fourier  transform  of  x.: 

K  J 


P .  =  l  X  exp(i2x1ik/N),  k  =  0,1 . N-l. 

*  j=0  -1 


(5-8) 


In  order  to  show  that  the  samples  p^  have  the  correct  corre¬ 
lation  properties,  the  autocorrelat ion  of  p^  can  be  written  as' 

N-l  N-l 

<Pk  P,*>  *  l  l  <x-  x  *>  exp[i2ir(jk-tm)/N]  .  (5-9) 

*  *  j=0  m=0  J  m 


However,  because  the  frequency  samples  Xj  are  independent, 


where 


<xj  V> a 


<C .  K  *>  =  2^ 


if  ,i  *  m 


if  j  *  m 


(5-10) 


(5-11) 


Therefore 


<pk  pt*>  =  2o2  |o  Sj  exPri2»j(k-t)/rt]  «  2o2  p(k-t)  (5-12) 

where  p(k-t)  is  the  discrete  autocorrelation  function  which  is  equal  to 
the  discrete  Fourier  transform  of  Sj. 
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The  first  100  decorrelation  times  of  the  received  or  combined 
power  for  the  one-way,  monostatic,  bistatic,  and  dual  channel  summation 
algorithm,  and  dual  channel  maximum  algorithm  cases  are  plotted  in  Figures 
13  a,b,c,d,e  respectively  for  sequences  generated  with  10  amplitude 
samples  per  one-way  decorrelat ion  time.  From  the  figures,  it  can  be  seen 
that  during  the  time  plotted  the  one-way  received  power  has  one  fade  more 
than  -25  dB  below  the  mean  power  and  none  more  than  -30  dB  below  the  mean 
power;  the  monostatic  received  power  has  ten  fades  more  than  -30  dll  below 
the  mean  power;  and  the  bistatic  received  power  has  two  fades  more  than 
-30  dB  below  the  mean  power.  The  dual  channel  combined  power  sequences 
have  no  faaes  more  then  20  dB  below  the  mean  power.  The  summation  algo¬ 
rithm  sequence  has  seven  fades  which  are  more  than  10  dB  below  the  mean 
power  but  none  of  those  fades  reaches  15  dB  below  the  mean.  The  maximum 
algorithm  has  eight  fades  which  fall  more  than  10  dB  below  the  mean  power 
and  two  that  exceed  15  dB  below  the  mean  power. 


Figure  13a.  Example  of  the  received  power  for  the  one-way  case. 
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Figure  13e.  Example  of  the  combined  power  for  identically  distributed 
dual  channels  (maximum  algorithm). 


/ 


The  phase  of  the  one-way  received  voltaqe. 


6(t)  =  tan-1{ Im[p(t)]/Re[p(t)]} 


(5-13) 


is  shown  in  Fiqure  14  where  it  can  be  observed  that  the  most  rapid  chanqes 
in  the  phase  occur  when  the  received  power  is  in  a  deep  fade.  The  auto¬ 
correlation  of  the  one-way  voltaqe  shown  in  Fiqure  13a  is  plotted  in  Fiq¬ 
ure  15  alonq  with  the  Gaussian  form.  Althouah  the  autocorrelation  of  one 
sequence  is  a  random  variate  and  therefore  has  a  distribution  about  its 
mean,  the  averaqe  autocorrelation  of  twenty  sequences,  plotted  in  Fiqure 
16,  shows  close  aqreement  with  the  desired  Gaussian  form. 
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Figure  14.  Phase  of  the  voltage  for  the  one-way  received  power  shown  in 
Figure  13a. 
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Figure  15.  Autocorrelation  of  one  random  sequence  of  the  one-way  received 
voltage  and  the  desired  Gaussian  form. 


Figure  16.  Average  autocorrelation  of  20  random  sequences  of  the  one-way 
received  voltage  and  the  desired  Gaussian  form. 


The  statistics  of  the  duration  and  separation  of  fades  or  flares 
is  gathered  by  stepping  through  the  random  sequences  and  comparing  the 
received  power  with  the  desired  level.  Because  of  the  finite  length  of 
the  voltage  sequences,  some  boundary  effects  exist.  For  example,  when  a 
new  sequence  is  generated  the  first  fade  (flare)  below  (above)  some  level 
L  can  be  recorded  only  after  encountering  the  first  amplitude  above 
(below)  L.  Similarly,  a  fade  or  flare  must  end  before  the  end  of  the 
sequence  to  be  included  in  the  data.  Also,  except  when  the  sequence  ends 
during  a  fade  or  flare,  the  number  of  intervals  between  fades  or  flares  in 
a  sequence  will  be  one  less  than  the  number  of  occurrences  of  fades  or 
flares.  Histograms  are  generated  of  the  samples  of  duration  anji  separa¬ 
tion  from  which  median  values  are  calculated.  In  addition,  statistics  are 

collected  on  the  amplitudes  z(t)  and  on  the  number  of  crossings  of  the 

i 

level  L. 

The  total  number  of  received  voltage  samples  generated,  which  is 
equal  to  the  number  of  voltages  per  sequence  times  the  number  of  sequences 
generated,  is  selected  in  order  to  achieve  about  104  fades  or  fjares.  The 
mean  number  of  occurrences  per  sequence  is  <N(L,t1 )>/2(N/n)  where  N/n  is 
the  number  of  decorrelation  times  per  sequence.  The  number  of  Samples  n 
per  decorrelation  time  has  been  selected  so  that  the  mean  duration  of  a 
fade  or  flare  is  roughly  greater  than  10  samples  in  order  to  collect  mean¬ 
ingful  duration  data.  Thus  for  very  deep  fades  with  correspondingly  short 
durations,  20  or  more  samples  per  decorrelation  time  are  required  to 
obtain  good  statistics  on  the  fade  duration.  Consequently,  the  number  of 
decorrelation  times  per  sequence  is  reduced  and  the  required  number  of 
sequences  that  must  be  generated  is  increased.  Because  of  these  consider¬ 
ations,  the  minimum  level  relative  to  the  mean  power  for  fades  has  been 
limited  to  -15  dB  and  the  maximum  level  for  flares  has  been  limited  to  +5 
dB. 
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SECTION  6 

SIMULATION  RESULTS 


The  results  from  the  Monte  Carlo  simulation  with  a  Gaussian  PSD  \  \ 

are  presented  in  this  section  for  the  one-way,  monostatic,  and  bistatic 
cases,  and  for  the  dual  channel  communication  link  cases.  First,  the 
moments  of  the  amplitude  z  for  each  of  the  cases  are  compared  with  the 
expected  values  given  in  Table  I  (Page  24)  in  order  to  verify  that  the 
simulated  received  voltages  have  the  desired  amplitude  statistics.  Then 
the  correlation  properties  of  the  received  voltages  are  verified  by 
comparing  the  mean  number  of  level  crossings  and  the  mean  duration  and 
mean  separation  of  fades  and  flares  with  the  results  calculated  in  Section 
4.  Finally,  the  mean,  median,  and  standard  deviation  of  the  duration  and 
interval  between  fades  or  flares  are  presented  in  tabular  form.  A  compar¬ 
ison  is  made  with  multiple  phase  screen  calculation  results  in  Appendix 
C  Plots  of  the  probability  density  function  and  the  cumulative  probabil¬ 
ity  distribution  are  shown  for  the  one-way  case  in  Appendix  D. 

The  first  four  moments  of  the  amplitude  z  of  the  received  volt¬ 
age  normalized  to  the  expected  values  are  snown  in  Table  II  for  the 
various  propagation  cases  qiven  in  Table  I.  The  largest  deviation  of  the 
amplitude  moments  from  the  expected  values  occur  in  the  monostatic  and 

/ 

bistatic  cases  while  the  one-way  amplitude  moments  vary  by  at  most  4.3  ./ 

percent  from  the  expected  values. 

The  statistics  of  one-way,  monostatic,  and  bistatic  fades  and 
flares  are  summarized  in  Tables  III,  IV,  and  V  respectively.  The  mean. 


56 


median,  and  standard  deviation,  normalized  to  t0  of  fade  or  flare  dura¬ 
tion  and  separation  are  shown  for  various  levels  along  with  the  mean-to- 
standard  deviation  (u/o)  ratios. 


Table  II.  Moments  of  the  received  voltages  from  the  Monte  Carlo 
simulation. 


Case 

<z">/t|7"]* 

No.  of  Amplitude 

n  ■  1 

n  •  2 

n«3 

n  ■  4 

One-Way 

1.007 

1.017 

1.079 

1.043 

344064 

Nonost at 1c 

1.010 

1.078 

1.049 

1.071 

574788 

Blstitlc 

0.995 

0.989 

0.935 

0.983 

574788 

Ouil  Channel  Coamunlcation 

Sumnitlon  Algorithm 

1.000 

1.001 

1.00? 

1.003 

574788 

Maximum  Algorithm 

0.997 

0.995 

0.993 

0.99? 

574788 

*<?">  denotes  simulation  results  for  0  dB  fading. 
E[  ?°]  denotes  expected  values  given  in  Table  I. 


Close  agreement  is  seen  between  the  analytic  results  of  Section 
4  and  the  simulation  results  for  the  mean  number  of  level  crossings  and 
the  mean  duration  and  separation  of  fades  and  flares.  The  close  agreement 
of  the  analytic  and  simulation  mean  values  confirms  that  the  simulated 
received  voltage  has  both  the  proper  amplitude  statistics  and  correlation 
properties  because  the  mean  number  of  level  crossings  depends  depends  on 
both  the  first  and  second  order  statistics  of  the  received  voltage. 

The  mean-to-standard  deviation  ratios  of  the  data  range  from  1.1 
to  1.8  with  most  in  the  range  from  1.2  to  1.5.  As  a  point  of  comparison, 
a  Rayleigh  distribution  has 

v/o  =  (4/it  -  l)-1/2  *  1.913  .  (6-1) 
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Table  V:  Statistics  of  bistatic  scintillation. 


Thus  the  durations  and  separation  have  narrower  distributions  about  their 
means  than  does  a  Rayleigh  distribution. 

The  statistics  of  dual  channel  fades  are  summarized  in  Table 
VI.  In  this  case,  no  analytic  expressions  are  available  for  the  mean 
quantities.  When  comparing  the  summation  algorithm  (denoted  by  SUM  in  the 
table)  with  the  maximum  algorithm  (denoted  by  MAX),  the  mean  durations  and 
separations  of  fades  are  slightly  larger  for  the  summation  algorithm  at 
all  levels  shown  while  the  mean  number  of  level  crossings  is  slightly 
smaller.  In  comparison  with  the  single  channel  one-way  communication  link 
data,  the  dual  channel  has  on  the  average  shorter  fade  durations  which  for 
deep  fades  below  -10  dB  occur  much  farther  apart.  Fades  -10  dB  below  the 
mean  power  in  the  dual  channel  system  occur  only  about  1/3  as  often  as 
they  occur  in  the  single  channel  system  and  they  are  only  2/3  as  long. 


Table  VI:  Statistics  of  dual  channel  fades. 


s 

u 

* 

H 

a 

X 

1/<S> 

<•»») 

n* 

<8(1. .n)> 

|  Duration  of  fades 

|  Separation  of  Fades  I 

S  taulat  Ion 

ho* 

Hr  an* 

Std.  Dev.* 

Had  I  an* 

n/o 

No* 

Hean* 

Std.  Dev.* 

Median* 

m/o 

X 

0 

10 

0.66? 

17339 

1.83 

1.60 

1.30 

1.14 

17311 

3.01 

1.78 

7.50 

1.69 

X 

-J 

TO 

0.674 

10599 

0.8?1 

0.663 

0.60 

1.74 

10546 

7.96 

1.94 

7.45 

1.53 

X 

-s 

» 

0.508 

9843 

0.556 

0.437 

0.433 

1.77 

971? 

3.90 

7.98 

3.03 

1.31 

X 

•10 

30 — 

0.143 

10?30 

0.771 

0.196 

0.70 

1.38 

9716 

13.3 

17.5 

9.43 

1.06 

X 

0 

!0 

0.633 

16581 

1.88 

1.57 

1.30 

1.70 

16559 

3.15 

1.81 

7.60 

1.74 

X 

-J 

TO 

0.6?? 

10436 

0.859 

0.653 

0.65 

1.3? 

10376 

3.70 

1.99 

7.65 

1.61 

X 

-s 

30 

0.46? 

9707 

0.590 

'  0.471 

0.467 

1.40 

9573 

4.28 

3.35 

3.73 

1.78 

X 

-10 

30 

0.l?9 

17896 

0.786 

0.194 

0.733 

1.47 

17175 

14.6 

13.9 

10.3 

1.05 

*Noraalt{ed  to  t0 

*Niaber  o f  samples  of  the  duntfon  or  separation 
**«uW>er  of  voltage  samples  per  t0 


SECTION  7 

RESULTS  ANd  CONCLUSIONS 


The  major  results  of  this  study  are: 

•  Binary  error  rates  under  slow  Rayleigh  fading  conditions  for 
various  frequency-shift  keying  and  phase-shift  keying  demodu- 

i 

lation  techniques  have  been  calculated  for  two-way  and  dual 
channel  communication  links. 

•  Analytic  expressions  have  been  derived  for  the  mean  number  of 
level  crossings  and  the  mean  duration  and  separation  of  fades 
and  flares  for  one-way,  monostatic,  and  bistatic  propagation 

j 

with  Rayleigh  amplitude  statistics  and  arbitrary  power  spec¬ 
tral  densities.  I 

! 

i 

•  The  probability  density  functions  and  cumulative  distribu¬ 
tions  of  the  duration  and  separation  of  fades  and  flares  for 
one-way  propagation,  Rayleigh  amplitude  statistics,  and  a 
Gaussian  power  spectral  density  are  compiled  in  Appendix  D. 

The  key  conclusions  of  this  study  are: 

•  Under  slow  Rayleigh  fading  conditions,  a  b’static  communica¬ 
tion  link  through  a  satellite  transponder  requires  about  9  dB 
more  bit  energy-to-noise  density  than  does  a  one-way  link  to 
achieve  a  bit  error  rate  of  10"3  while  a  dual  channel  commun¬ 
ication  link  requires  about  10  dB  less  bit  energy-to-noise 
density  to  achieve  a  10"3  bit  error  rate. 
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Relative  to  monostatic  space  based  radar  systems,  bistatic 
radars  have  shorter  mean  fade  durations  and  longer  mean 
separations  of  deep  fades  more  than  10  dB  below  the  mean 
power. 

Fades  in  dual  channel  communication  links  which  are  more  than 
10  dB  below  the  mean  power  are,  on  the  average,  two-thirds  as 
long  and  occur  about  one  third  as  often  as  fades  10  dB  below 
the  mean  power  in  single  channel  links. 
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APPENOIX  A 

SCALING  FACTORS  FOR  fv  POWER  SPECTRAL  DENSITIES 


It  was  shown  in  Section  4  that  the  mean  number  of  level  cross¬ 
ings  and  the  mean  duration  and  separation  of  fades  or  flares  for  an 
arbitrary  power  spectral  density  (PSD)  can  be  obtained  from  the  mean 
values  for  a  Gaussian  PSD  by  simple  scaling  factors.  It  is  the  purpose  of 
this  Appendix  to  calculate  the  scaling  factors  for  general  f"w  power 
spectral  densities  and  to  sketch  the  formalism  for  calculating  the  scaling 
factors  for  arbitrary  power  spectral  densities. 

In  order  to  make  these  calculations  as  simple  as  possible,  it  is 
assumed  in  bistatic  geometries  that  the  two  one-way  paths  have  identical 
first  and  second  order  statistics.  It  is  an  easy  matter  to  generalize 
these  results  to  a  bistatic  geometry  with  different  second  order  statist¬ 
ics  on  each  path  (i.e.,  to  geometries  where  the  PSD  and/or  the  decorrela¬ 
tion  time  of  the  two  paths  differ). 

A.l  GENERAL  RESULTS 

8efore  calculating  the  scaling  factors,  some  general  results 
need  to  be  established.  The  autocorrelation  function  of  the  complex  volt¬ 
age  p (t )  is  related  to  the  PSD  S(f)  by  the  Fourier  transform  pair 

p(t)  *  /  S(f)e^*  fTdx  ( A-l) 
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s(f)  =  I  p(t)e 

.00 


i2irfTdt 


(A-2) 


For  a  Gaussian  autocorrelation  function 

p( t )  =  exp[-(t/ tj ) 2  ]  (A-3) 

it  is  well  know  that  the  PSD  is  also  Gaussian  (c.f.  Equation  5-1). 

The  scale  factor  A,  defined  in  Equation  4-4,  is  proportional  to 
the  second  time  derivative  of  p(t)  evaluated  at  zero  time  delay. 
Differentiating  Equation  A-l  twice  qives 

p( o )  =  -4 n2  /  f2S(f)df  (A-A) 


2 

which  for  a  Gaussian  PSD  reduces  to  -Zh^.  The  scale  factor  A  for  an 
arbitrary  PSD  is  then 


A2  =  2i»2Tj2/  f2S(f )df  . 


( A- 5 ) 


Convergence  of  the  integral  in  Equation  A-5  requires  that  S(f)  approach 
zero  faster  than  f-3  for  large  f. 


The  other  scale  factor  5,  defined  in  Equation  2-26,  relates  the 
two-way  decorrelat ion  time  to  the  one-way  decorrelat ion  time.  It  was 
shown  in  general  in  Section  2  that  for  two-way  geometries 


p.  (  t) 

two-way 


p2  (  t) 

one-way 


(A-6) 


'  "  i 

•  •  • 

« 
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The  one-way  decorrelation  time  tj  is  defined  by 


p  (it  )  =  1/e  • 

one-way  1 


(A- 7) 


The  decorrelation  time  of  the  received  voltage  in  a  two-way  propagation 
geometry,  t0.  is  then  calculated  from  the  equation 


p  (t0)  =  l//e  . 

one-way  0 


(A-8) 


Again  for  a  Gaussian  PSD,  t0  is  just  x^Z/2  so  that  in  general 


6  =  /?T  /t  i  • 

o  1 


(A-S) 


A, 2 


is 


f-*  (u  >  3)  POWER  SPECTRAL  DENSITIES 

A  simple  form  for  the  PSD  with  an  f*u  frequency  dependence 


S(f) 


_ So _ 

(f2  +'fo2)M/2 


(A-10) 


where  the  coefficient  S0  is  given  by  the  normalization  condition 


p(o)  «  /  S(f )df  =  1 


( A- 1 1 ) 


which  reduces  to 

<  ..  fo  u'1  r(M  12) 

0  r (1/2)  r[ (p-l)/2]  • 
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( A— 12) 


••  i 
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Note  that  Equation  A-10  is  valid  only  for  p  >  3  because  of  the  requirement 
that  p(o)  be  finite.  A  second  more  general  form  with  an  inner  scale  which 
is  valid  for  p  3  is  considered  later  in  this  appendix.  The  autocorrela¬ 
tion  function  corresponding  to  the  PSD  of  Equation  A-10  is 


p(t) 


(2nf0T) 


(u-i  )h 


K 


(u-O/2 


(2*foT) 


2(w_3  )/z  r[{u-l)/2}] 


(A-13) 


where  ICy  is  the  vth  order  K  Bessel  function.  The  one-way  and  two-way 
decorrelation  times  can  now  be  calculated  by  solving  the  transcendental 
equations 


((u-l)/2  K 


(u-i  )/2 


(x) 


.  1 


2(m-s  )/2 


( A— 14 ) 


Jn-1  )/2 


K(m-x)/2  {y) 


.  1 


2(m-3  )/2  r[  (m-1)/2]  fe 


( A— 15 ) 


for  x  and  y  where 

x  =  2*  fQti 

( A— 16 ) 

y  =  ?1,foTo  • 

( A— 1 7) 

The  scale  factor  6  is  simply 


6  =  /2y/x 


(A-18) 
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which  is  only  a  function  of  y  .  The  scale  factor  A  is  given  by  evaluating 
Equations  A-5  and  using  Equation  A-16  to  write  f0  in  terms  of  x  and  tj . 
This  results  in 


a2  .  r[(u-3)/2l  x2 
4  r[ (p-l )/2] 


(A-19) 


which  is  again  only  a  function  of  y  and  is  therefore  independent  of  the 
decorrelation  time  and  of  the  level  of  the  fades  or  flares. 


The  scale  factor  5  is  plotted  in  Figure  A-l  for  3  <  y  <  10  where 
it  can  be  seen  that  o  approaches  unity  as  u  increases.  Recall  that  6  and 
A  are  both  unity  for  a  Gaussian  PSD.  The  scaling  for  the  mean  duration 
and  separation  of  fades  or  flares,  1/6A  (c.f.,  Equation  4-29),  is  shown  in 
Figure  A-2  for  one-way  (5  =  1)  and  two-way  propagation.  For  y  =■  4,  the 
scale  factors  are  as  follows: 


6  =  0.8946 


1/6A  = 


0.6589  one-way 
0.7365  two-way 


(A-20) 


Note  that  1/6A  changes  rapidly  between  f4  and  f”3  where  1/6A  is  zero. 


f-w  POWER  SPECTRA!.  DENSITIES  WITH  AN  INNER  SCALE 


The  PSD  of  Equation  A-10  can  be  used  for  y  >  3  only.  A  more 
general  form  has  been  given  by  Shkarofsky13  as 


S(f)  = 


<w1/z  y2  lay*2  ♦  H 

V-o/2  (2*fo'o>  *  f|i“/2 


(A-21) 


which  depends  on  both  an  outer  scale  f  and  an  inner  scale  t”.  This 

o  o 

PSD,  which  is  shown  schematically  in  Figure  A-3,  is  valid  for  any  y.  For 
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log(f) 


Figure  A-3.  Schematic  diagram  of  the  power  spectral  density  with  an  inner 
and  an  outer  scale. 

frequencies  in  the  range  fQ  <<  f  «  t"1  ,  S(f)  has  an  approximate  f"u 
frequency  dependence.  Once  y  has  been  fixed,  the  outer  scale  f0  is 
determined  by  the  decorrelation  time.  The  other  scale  is  left  as  a  param¬ 
eter.  Because  an  f“3  PSD  provides  a  slow  scattering  limit  to  Rayleigh 
fading  environments,  the  subsequent  discussion  will  consider  only  y  =  3. 
Setting  y  =  3  in  Equation  A-21,  the  corresponding  autocorrelation  function 
is 


where 


p(t) 


/b2  +  ( 2-it f qt )Z  Kj[/b2  +  ( 2ir f qt )2 ] 

Bki  (a) 


B  =  2*  ft  . 
0  0 


(A-22) 


(A-23) 
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In  this  case,  the  one-way  and  two-way  decorrelation  times  are  computed  by 
solving  the  equations 


42  +  x2  Ki  [/e2  +  x2 
SKX  (0 ) 


1 

e 


/b2  +  y2  [/b2  +  /]  _  _i 
bMb)  ~/i 


for  x  and  y  where 

x  =  2*  f  0t  i 


{ A— 24) 

(A-25) 


(A-26) 


y  *  2rtf0T0 

and  where  x  and  y  are  now  functions  of  the  parameter  b.  Once  x  and  y  have 
been  obtained.  Equation  A-18  is  used  to  determine  6  as  a  function  of  b. 

The  scale  factor  a  is  determined  by  evaluating  Equation  A-5  which  yields 


a2  =  Mb)x2 

6Mb) 


(A-28) 


The  scale  factor  6  and  the  scaling  for  the  mean  duration  of  fades  and 
flares  (1/5A)  are  plotted  in  Figures  A- 4  and  A-5  respectively  with  b  in 
the  range  lO"10  <_  6  <_  1  .  Note  however,  that  as  b  approaches  unity,  the 
range  of  frequencies  for  which  the  PSD  has  an  f-3  frequency  dependence 
shrinks  to  zero.  For  B  <  10"3 ,  Figure  A-5  shows  that  1/6A  is  only  weakly 
dependent  on  8,  and  for  b  <  10"2 ,  Figure  A-4  shows  that  5  is  independent 
of  b  . 


All  that  remains  to  complete  this  discussion  is  to  relate  B  to 
physical  quantities.  To  this  end.  Equation  A-26  is  used  to  related  f0 


74 


to  the  decorrelation  distance  and  the  solution  of  Equation  A-24.  The 
inner  scale  t*1  can  be  estimated  by  noting  that  physically  the  PSD  must  go 
rapidly  zero  for  frequencies  that  exceed  the  radio  transmission  frequency 
fpp.  The  inner  scale  is  then  n/fpp  where  n  is  a  parameter  which  is 
the  order  of  unity.  The  value  of  3  is  then  determined  by 

e/x(e)  =  n/ ( t  ^  f  pp  5  (A-29) 

where  x( e )  is  the  solution  of  Equation  A-24  which  is  itself  a  function  of 
B.  This  equation  is  easily  inverted  by  plotting  &  versus  x(b)/b  which  is 
equal  to  TjfRp/n.  The  results  are  shown  on  Figure  A-6.  Because  x(b) 
for  b  <  10'2  is  independent  of  3  and  constant  at  a  value  of  1.658,  B  is 
accurately  given  by 

B  =  1.658n/(TlfRF),  B  <  10‘2  .  (A-30) 


Figure  A-6.  0  as  a  function  of  and  the  transmission  frequency  f^p. 
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The  parameter  g  is  a  function  of  the  decorrelation  time  and  hence  the 
scaling  factors  are  also  functions  of  tj  . 

The  scaling  factors  for  satellite  systems  with  transmission 
frequencies  in  the  range  250  KHz  to  23  GHz  are  shown  in  Table  VII  for  slow 
fading  conditions.  Typically,  the  one-way  scaling  factor  1/a  is  0.18  and 
the  two-way  scaling  factor  1/6A  is  0.21. 


Table  VII.  Slow  fading  (f“3  PSD)  scaling  factors  for  satellite  systems. 


fRF(C-Hz) 

i!  (sec) 

log{0  )* 

6 

1/6A 

one-way 

two-way 

0.25 

20 

-9.5 

0.844 

0.182 

0.216 

1.2 

10 

-9.9 

0.844 

0.178 

0.211 

7.5 

3 

-10.1 

0.844 

0.176 

0.209 

20.0 

2 

-10.4  I  0.844 

0.174 

0.206 

-10.4 


0.844 


0.174 


0.206 
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APPENDIX  B 

LEVEL  CROSSING  PROBLEM  FOR  BISTATIC  PROPAGATION  WITH 
INDEPENDENT  SECOND  ORDER  STATISTICS 


The  mean  number  of  level  crossinqs  is  calculated  in  this  Appen¬ 
dix  for  bistatic  propaaation  where  the  second  order  statistics  of  the  two 
paths  are  independent.  The  amplitudes  of  the  one-way  voltaqes  of  the  two 
paths  are  assumed  to  be  Rayleiqh  distributed  and  independent. 

Let  0 j ( t )  and  p2(t)  denote  the  autocorrelation  functions  of  the 
two  paths.  Both  the  power  spectral  densities,  or  equivalently  the  func¬ 
tional  form  of  the  autocorre 1  at  ion  functions,  and  the  decorrelation  times 
are  allowed  to  differ.  The  calculation  of  the  mean  number  of  level  cross¬ 
ings  then  follows  the  deviation  of  Section  4. 

The  joint  probability  density  function  of  the  amplitude  and  the 
time  derivative  of  the  amplitude  for  the  i**1  one-way  propagation  path  is 
given  by  Equation  4-11  as 


r  exp(-r2/?o2)  expf-ij jr z/4A^a2) 


( i  *  1  or  2) 


where  is  defined  in  Equation  4-4  and  r l ^ is  the  one-way  decorrela 
tion  time  of  the  i*h  path.  The  joint  probability  density  function  of 
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the  amplitudes  s  and  t  of  the  two  independent  paths  and  of  the  time 
derivatives  of  s  and  t  can  be  written  as 


!\  V 


O14  (4hA1A20Z/t1  ^  tTi  >2) 

(8-2) 


Usinq  the  same  chanqe  of  variables  that  was  used  in  Section  4,  the  joint 
probability  density  function  of  the  received  amplitude  z  and  the  time 
derivative  z  becomes 


:u 


}  . 
*  r 


f  ( 2 ,  Z )  «  (t1i1T1i2/4-.A1A206)  /  ds  /  ds  2S- 


(B-3) 


exp  | 


.[s2>z2/s2+tJi1sz/2a5+t5>2(z2s-2+zzs2s-u-?zzss'3)/2a2]/Zo2}  . 


Only  the  s  inteqral  can  be  performed  in  closed  form  which  reduces  Equation 
B-3  to 


where 


f(z,z)  =  (t^z/2/kAjo5)  /  ds  {csu^z2)-1/2 


x  exp i  -[s2+z2s‘2+tf^ zcs2(?sl,+z2  )-1z2/2Aj  ]/2o2} 


(8-4) 


'•>7 


'  I 


\f/ 


'!■  J ' 

v-  ?  • 

;  v.;;"  / 


I  ;»  *) 


-  i 


:  ./■ 
i/.7' 


2  2  2  2 
C  =  *1  iV(Tl  2A;) 
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(B-5) 


v 

/ 


-  r\f:  \ 


i  a-  ( 


.i  •  v' 


\  /  . 


■  n  /  . 

n/  v 

>  . 

\  f ;  ■ 
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The  mean  number  of  level  crossings  is  obtained  by  substituting  Equation 
B-4  into  Equation  4-3  and  performing  the  integral  over  z.  After  changing 
the  variables  in  the  remaining  integral  to  u  =  s/a,  using  Equation  2-12  to 
write  o'*  in  terms  of  the  mean  received  power  <S>,  and  using  Equation  2-26 
to  write  tj  in  terms  of  the  decorrelation  time  of  the  received  signal, 
the  mean  number  of  crossing  of  the  power  level  L  becomes 


<N(l,T)>  =  6ax(J/,o) 


-|1  /2 

16L/if  <S>  ' 

1+(T1  ,l/Ti  ,2  )2. 


«  /  /c+4L/<S>u“  exp{  -[  l+4L/<S>u1*]  u2/2}  du 

o 


( B-6) 


This  expression  reduces  to  Equation  *-26  when  Tl>1  =  |t1j2  and  Aj  =  a2  . 


i 

1 
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APPENDIX  C 
MPS  RESULTS 


The  first  four  moments  of  the  amplitudes  of  the  five  KPS  reali¬ 
zations  that  have  been  analyzed  for  this  report  are  presented  in  Table 
VIII.  Each  realization  consists  of  16384  amplitude  samples  and  the  corre¬ 
sponding  phases.  The  data  have  been  normalized  so  that  the  average  power  / 

is  unity.  Hence  <z2>  is  identically  one  in  each  case.  The  other  moments 
show  close  agreement  with  Rayleigh  values. 


*<zn>  denotes  the  average  from  16384  MPS  samples 
E[  zn]  denotes  the  Rayleigh  value 


83 


t  • 


The  realizations  differ  in  the  number  of  samples,  n,  per  decor¬ 
relation  time  which  varies  from  72.8  for  MPS  8091  to  3.35  for  MPS  8143. 

In  the  latter  case,  the  statistics  of  deep  fade  durations  will  be  unreli¬ 
able  because  almost  all  of  the  deep  fades  will  last  for  just  one  sample. 
With  this  in  mind,  the  means,  medians,  standard  deviations,  and  mean-to- 
standard  deviation  ratios  are  presented  in  Tables  IX  through  XIII  for  the 
five  MPS  rea  1  i  z  at  ions . 

Table  IX:  Statistics  of  fades  for  the  MPS  8091  realization. 


Our At  ion  of  f ados 

No*  Mean*  Std.  |)pv.#  Median* 
11?  1.70  1.11  0.7/? 

117  P.64'»  0.M3  0  46 1 

iro  0.406  0.406  0.1S8 

81  O.?40  0.761  0.1/8 

44  0.1S1  0.118  0.006 

Nnrmj  1  i /od  to  i  0 

of  samples  of  the  durat  ton  or  separat  ton 
**Nun#'nr  of  voltage  samples  per  r0 


Table  X:  Statistics  rf  fades  for  the  MPS  8094  realization. 


^Normal 1/Pd  to  r0 

+Nutnher  of  samples  of  the  duration  or  separation 
M Number  of  voltage  samples  per  t0 
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Table  XI: 


Statistics  of 


fades  for  the  MPS  8143  realization. 


1/<S> 

<dB) 

Duration  of  Fades 

...  _ .  ___  . .  _ ...  _ : 

Separat ion  of 

Fades 

n»* 

<N(l.,0)> 

Nof 

. 

Mean4 

Std.  Dev.* 

Median* 

v/o 

tk>* 

Mean* 

Std.  Dev.* 

Median* 

m/o 

0 

3.35 

0.631 

1543 

1.99 

1.3? 

1.19 

1.16 

154? 

3.17 

1.09 

?.39 

1.67 

-3 

3.35 

o.;?5 

1335 

1.09 

0.034 

0.590 

1 .75 

1374 

?.76 

1.54 

?.09 

1.79 

-5 

3.35 

0.6R8 

1606 

0.00? 

0.603 

0.?99 

1.3? 

1605 

?.90 

1.00 

?.09 

1.61 

-10 

3.35 

0. 

1051 

0.460 

0.?90 

0.?99 

1.50 

1050 

4.66 

3.95 

3.79 

1.10 

-15 

3.35 

0.104 

45? 

0.355 

0.146 

0.799 

7.44 

451 

10.7 

i  . . 

10.3 

7.15 

1.04 

Normal  lied  to  t0 

Ntaaber  of  samples  of  the  duration  or  separation 
4Nim*er  of  voltage  samples  per  t0 


Table  XII:  Statistics  of  fades  for  the  MPS  8642  realization. 


1/<S> 

(dB) 

n‘* 

Dui  at  ion  of  Fades 

Separation  of 

Fades 

No4 

Mean4 

Std.  Dev,* 

Median* 

v/o 

Wo4 

Mean4 

.. 

Std.  Dev.* 

Median* 

ll/<l 

n 

3.40 

0.565 

1370 

7.74 

1.93 

1.44 

1.16 

137? 

3.54 

7.14 

7.87 

1.65 

-3 

3.40 

0.650 

1579 

1.71 

0.900 

0.574 

1.71 

1578 

3.08 

1.69 

7.30 

1.8? 

-5 

3.40 

0.615 

1447 

0.805 

0.69ft 

0.787 

1.77 

1446 

3.75 

7.0? 

7  58 

1.61 

-10 

3.40 

0.399 

939 

0.497 

0.347 

0.787 

1.43 

938 

5.01 

4.70 

3.44 

1.19 

-15 

3.40 

0.17? 

404 

0.371 

0.709 

0.707 

1.78 

403 

11.6 

17.0 

7.18 

0.967 

Norma  Hied  to  t0 

Number  of  samples  of  the  duration  or  separation 
4Numher  of  voltage  samples  per  t0 
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Table  XIII:  Statistics  of  fades  for  the  MPS  9069  realization. 


l/<$> 

(dB) 

Ouratton  of  fades 

Separation  of 

Fades 

n** 

No* 

Mean* 

Std.  Dev.* 

Med  1 4n* 

y/o 

Ho* 

Mean* 

Std.  Dev.* 

Hed i an* 

tt/o 

mm 

10.5 

0.554 

433 

?.30 

?.44 

1.34 

0.943 

432 

3.61 

2.60 

2.77 

1.39 

n 

10.5 

0.664 

514 

1.71 

1.08 

0.763 

1.1? 

518 

3.01 

1.93 

2.48 

1.56 

-5 

10.5 

0.616 

481 

0.9?3 

0.801 

0.572 

1.15 

480 

3.24 

2.28 

2.48 

1.42 

-10 

10.5 

0.44? 

3*5 

0.453 

0.354 

0.286 

1.28 

344 

4.5? 

4.60 

2.86 

0.983 

-15 

10.5 

0.?61 

?04 

0.259 

0.184 

0.095 

1.41 

203 

7.67 

8.93 

4.86 

0.859 

‘Normal  I  red  to  t0 

*Nu»4>er  of  samples  of  the  duration  or  separation 
**Nim4>er  of  voltaqe  samples  par  t0 


The  average  and  the  standard  deviation  of  the  MPS  scaling 
factors  for  the  mean  duration  and  separation  of  fades  and  flares  have  been 
calculated  from  the  ratio  of  the  mean  number  of  level  crossings  for  the 
MPS  realizations  to  the  mean  number  of  level  crossings  for  a  Gaussian 
power  spectral  density.  The  mean  scaling  factor  is  calculated  as 


<1/a  > 


i  j,  <nVVV<"V'-o>>mps 


(C-1) 


where  ,  i  =  1  through  5,  are  the  five  levels  of  the  tabulated  MPS 
data.  The  mean  duration  and  separation  data  are  not  included  in  the 
average  because  these  values  are  more  sensitive  to  the  number  of  amplitude 
samples  per  decorrelation  time  and  are  not  independent  of  the  mean  number 
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of  level  crossings.  The  results  are  shown  in  Table  XIV.  Realizations 
8143,  8642,  and  9069  all  have  average  values  of  1/a  of  about  unity  which 
indicates  that  their  power  spectral  densities  are  nearly  Gaussian.  The 
power  spectral  densities  of  the  other  two  realizations  can  be  estimated  by 
using  <1/a>  and  the  one-way  curve  of  Figure  A-2  to  find  an  effective 
f-P  dependence.  The  power  spectral  density  is  then  approximately 
f-3-8  for  MPS  8091  and  is  f"4 -7  for  MPS  8094. 


Table  XIV:  Average  Scaling  Factors  (1/a). 


Realization 

<1/A> 

Std.  Dev.  (1/A) 

8091 

WEM 

0.054 

8094 

wmm 

0.030 

8143 

0.240 

8642 

1.180 

0.238 

9069 

1.048 

0.015 

I 

( 
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APPENDIX  D 


PROBABILITY  DENSITY  FUNCTIONS  AND  CUMULATIVE 
PROBABILITY  DISTRIBUTIONS  FOR  ONE-HAY  SCINTILLATION 


The  probability  density  function  histograms  and  cumulative 
distributions  of  the  duration  and  separation  of  fades  or  flares  under 
Rayleigh  fading  conditions  with  a  Gaussian  power  spectral  density  are 
presented  here  for  oneway  propagation.  All  levels  are  relative  to  the 
mean  power.  In  a  few  cases  noted  on  the  figures,  the  data  have  been 
rebinned  to  give  the  probability  density  histograms  a  more  continuous 
appearance.  Otherwise,  each  interval  of  the  histogram  corresponds  to  one 
sample  of  the  received  voltage  sequences.  The  mean  (m),  median  (M),  stan¬ 
dard  deviation  (a),  (all  normalized  to  t0)  and  number  of  samples  are 
indicated  on  the  figures.  For  fades  15  dB  below  the  mean  power  and  flares 
5  d3  above  the  mean  power,  the  conflicting  conditions  of  short  durations, 
which  requires  a  large  number  of  samples  per  decorrelation  time,  and  long 
intervals  between  events  resulted  in  about  10  and  20  percent  of  the  separ¬ 
ation  data  respectively  exceeding  the  separation  histogram  array  size. 
However,  because  the  moments  were  calculated  from  the  raw  data  and  less 
than  20  percent  of  the  data  were  effected,  the  mean,  standard  deviation, 
and  median  values  are  uneffected  by  the  improper  array  sizes. 
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Figure  b~2.  Cumulative  probability  distribution  of  the  deration  of  fades  below  0  dB. 


Figure  0-5.  Probability  density  function  of  tbe  duration  of  fades  below  -3  d8. 
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lative  probability  distribution  of  the  duration  of  fades  below  -3  d3. 


Cumulative  probability  distribution  of  the  separation 


Probabi  1 ity  dens  i 


Cumulative  probability  distribution  of  the  duration  of  fades  below  -5  dB. 


Figure  D-12.  Cumulative  probability  distribution  of  the  separation  of  fades  below  -5  dP. 


the  duration  of  the  fades  below  >10  dB 


Figure  D-14.  Cumulative  probability  distribution  of  the  duration  of  fades  below  -10  dB. 
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lity  density  function  of  the  separation  of  fades  below  -10  dB. 


Figure  D-16.  Cumulative  probability  distribution  of  the  separation  of  fades  below  -10  dB. 
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Figure  0-17.  Probability  density  function  of  the  duration  of  fades  below  -15  dfl. 
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Cumulative  probability  distribution  oi  the  duration  of  fades  below  -15  dB. 


Figure  0-20.  Cumulative  probability  distribution  of  the  separation  of  fades  below  -15  d3 


Figure  D-22.  emulative  probability  distribution  of  the  duration  of  flares  &ove  0  dB. 
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Figure  0-23.  Probability 


Cumulative  probability  distribution  of  the  separation 


Figure  D-30.  emulative  probability  distribution  of  the  duration  of  flares  above  5  dB. 
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Figure  0-31.  Probability 
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